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Abstract 

The quantum Zakharov system is described in terms of a La- 
grangian formalism. A time-dependent Gaussian trial function ap- 
proach for the envelope electric field and the low-frequency part of 
the density fluctuation leads to a coupled, nonlinear system of ordi- 
nary differential equations. In the semiclassic case, linear stability 
analysis of this dynamical system shows a destabilizing role played by 
quantum effects. Arbitrary value of the quantum effects are also con- 
sidered, yielding the ultimate destruction of the localized, Gaussian 
trial solution. Numerical simulations are shown both for the semiclas- 
sic and the full quantum cases. 



1 Introduction 

Quantum plasmas have received much attention in recent times, especially 
because of the ongoing miniaturization of ultrasmall electronic devices and 
micromechanical systems pQ and to the relevance of quantum effects for in- 
tense laser-plasmas [2 J and for dense astrophysical objects [3|. Frequently, 
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the de Broglie wavelength of the charge carriers (electrons, positrons, holes) 
of these systems is comparable to the characteristic dimensions of the sys- 
tem, making a quantum treatment unavoidable. Advances in the area in- 
cludes construction of quantum ion-acoustic waves [I], quantum magnetohy- 
drodynamics theories [5], quantum beam instabilities [6] [8] and shear Alfven 
modes in ultra-cold quantum magnetoplasmas [§]. New quantum collective 
excitations have also been identified for ultra-cold dusty plasmas [TU]-[To]. 
where quantum effects can be used for plasma diagnostics. Recently, spin 
effects have been included for non relativistic quantum plasmas [HH [T7] . 
Possible decisive applications of spin effects in quantum plasmas can appear 
in solid state plasmas as well as in pulsars and magnetars, with very intense 
magnetic fields (greater than 10 8 T). A more detailed review on quantum 
plasma models and their range of validity can be found in |18j . 

The quantum Zakharov equations [19], the subject of the present work, 
form a set of coupled, nonlinear partial differential equations for the envelope 
electric field and the slow part of the density fluctuation in an electron-ion 
electrostatic quantum plasma. It models the interaction between quantum 
Langmuir and quantum ion-acoustic waves. Exactly as for the classical Za- 
kharov system [20], the derivation of the quantum Zakharov equations comes 
from the existence of two time scales, a fast one associated to quantum 
Langmuir waves, and a slow one associated to quantum ion-acoustic waves. 
Sample applications can be found for quantum decay and four-wave insta- 
bilities, with relevant changes of the classical dispersions [19j. The quantum 
Zakharov system was also analyzed for the enhancement of modulational in- 
stabilities due to combination of partial coherence and quantum corrections 
|21j . More recently, the coupling between nonlinear Langmuir waves and 
electron holes in Wigner-Poisson quantum plasmas was studied via a two 
time scales formalism [22J. 

The existence of coherent structures, as soliton solutions for instance, is 
a relevant issue for any system of evolution equations. As an example from 
quantum plasmas, stable vortices and dark solitons have been constructed 
for Schrodinger-Poisson quantum electron plasmas [23]. At quantum scales, 
the transport of information in ultracold micromechanical systems can be 
addressed by means of such nonlinear structures. The basic objective of the 
present work is the investigation of the quantum effects for the existence of 
localized solutions for the quantum Zakharov system. Unlike the approach 
of Yang et ai, where exact bright solitons, gray solitons, W-solitons and M- 
solitons were found for the quantum Zakharov system [21], here approximate 
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solutions are obtained through a variational formulation and a trial function 
method. Exact solutions are of course very relevant, but variational solutions 
provides more insight on the role of quantum effects. For instance, the clas- 
sical Zakharov equations admit the Langmuir soliton solution [25]. Using a 
Gaussian ansatz as a trial function extremizing an action functional, one can 
get information about the perturbation of the Langmuir soliton by quantum 
effects. A priori, one can expect that wave-packet spreading and tunneling 
tends to enlarge the width of localized wave solutions. Other possibility is 
the appearance of instabilities of pure quantum nature, eventually destroying 
any coherent structure. Besides these considerations, the construction of a 
variational formulation for the quantum Zakharov equations is important by 
itself. Notice that the internal vibrations of solitary waves for the classical 
Zakharov system were analyzed by a variational approach using Gaussian 
trial functions (26]. The present contribution extends this work to the quan- 
tum realm. Similar time-dependent variational methods were also used, for 
instance, for the nonlinear pulse propagation in optical fibers (27] and for 
Bose-Einstein condensates [28J. 

Variational methods can indicate a general tendency of s system for which 
no general closed form solution is available. For instance, one can study 
the changes of localized or solitonic trial functions under the changes of a 
control parameter. The quantum Zakharov equations (see Section II) possess 
a single dimensionless quantity H measuring the importance of quantum 
perturbations. It is one basic task of this work, to analyze the changes in 
Gaussian trial function solutions for the quantum Zakharov system induced 
by modifications in H. 

This work is organized as follows. In section II, the quantum Zakharov 
system is described by a variational formulation. A variational solution in 
the form of a Gaussian ansatz is then proposed, in order to reproduce the 
main properties of the Langmuir soliton solution admitted in the classical 
limit. This time-dependent trial function approach leads to a dynamical sys- 
tem which can be analyzed for several parameter regimes. In section III, only 
the first-order quantum correction is retained, yielding a set of two coupled, 
nonlinear second-order ordinary differential equations for the widths of the 
envelope electric field and density perturbation. This nonlinear system is an- 
alyzed for its linear stability properties as well as for the existence of bounded 
solutions. In section IV, arbitrary strength of the quantum effects is allowed, 
resulting in a full system of equations. Further, the ultra-quantum case where 
quantum effects are the more relevant influence is analyzed, showing the ulti- 
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mate destruction of the Langmuir soliton due to wave-packet spreading and 
tunneling. Section V is reserved to the conclusions. 



2 Variational formulation 



The one- dimensional quantum Zakharov equations reads [19 
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where E = E(x, t) is the envelope electric field and n = n(x, t) is the density 
fluctuation. All quantities are expressed in a convenient dimensionless form. 
Further, 

HuJa 

K>B J-e 

is a parameter expressing the ratio between the ion plasmon energy and the 
electron thermal energy, where % is the scaled Planck constant, Kb the Boltz- 
mann constant, Ui the ion plasma frequency and T e the electron temperature. 
The formal classical limit is obtained for H = 0, yielding the original Za- 
kharov system. For more details on the derivation of the system ([IH2]) as 
well as for sample applications, see [T9"] . 

The one- dimensional quantum Zakharov equations are derived from the 
Lagrangian density 
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where it was introduced the auxiliary variable u so that 

du 
dx 

Indeed, 
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The last equation reproduces ([2]) after differentiation with respect to x. 

The classical Zakharov system is not integrable. However, it admits [25] 
the exact Langmuir soliton solution 

E = Eosech^exp^) , (9) 

n = -£ 2 sech 2 (^) , (10) 

where E is an arbitrary real parameter. Strictly, collisions of Langmuir "soli- 
tons" does not simply imply phase shifts between them [25], as expected for 
solitonic objects. Nevertheless, it is the interplay between nonlinear and dis- 
persive terms in the classical Zakharov system which allows for the existence 
of such coherent structures. Physically, the Langmuir soliton represents a 
hole in the low frequency part of the electron-ion density maintained self- 
consistently by the ponderomotive force. It would be desirable to achieve a 
better understanding of the quantum effects for this soliton solution. A priori 
one can expect tunneling of electrons trapped in the self-consistent potential, 
perturbing the remarkable stability of the Langmuir solitons in the classical 
case. Indeed, isolated classical Langmuir solitons do not decay [25] . 

A variational solution which reproduces the gross features of fl9l fT0|) and 
is at the same time analytically accessible is the Gaussian ansatz 

E = Aexp(-— + i(p + iKx 2 ), (11) 

n = -5exp(-^), (12) 

where A, B, a, b, cf) and k are functions of time. We assume A and B positive 
to maintain resemblance with (l9l fT0l) . As for any time- dependent variational 
method, notice that a main drawback of the Gaussian ansatz is that it does 
not allow for changes in the shape of the solution. The classical Zakharov 
system can be treated by a variational approach using a combination of 
Jacobi elliptic functions [29], but we use Gaussian functions for the sake of 
simplicity. 
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In order to calculate the Lagrangian L = ff^ Cdx corresponding to (TTTT - 
[T2j) there is the need of the derivatives du/dx and du/dt of the auxiliary 
function u. Combining (jSJ) and ffl2l) and introducing M = Bb, it follows that 



du 

at 
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b 2 



exp( 



x 
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T 



£eM- y ¥ )dy. 



(13) 



Inserting the last expression into the integral for the Lagrangian one con- 
cludes that it converges if and only if M = 0, implying a restriction on the 
allowable variational density functions. Indeed, if M ^ 0, then the proposed 
Gaussian ansatz leads to divergence of the Lagrangian. 

M being invariant is consistent with the conservation of the low frequency 
part of the mass 

1 



. ndx = Bb, (14) 
yix J-oo 

the last equality following from the variational solution and with the factor 
l/v^r being introduced for convenience. In addition, we note the existence 
of the conservation of the number of high frequency quanta 
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The numerical factor I/^/tt was used for convenience again, in the definitions 
of N and H, while the last equalities at (115H16I) follows from the proposed 
variational solution. The quantum Zakharov equations also preserve a mo- 
mentum functional, but this information is useless in the remaining. 



6 



Now we get the effective Lagrangian 

r- r ,2 l A 2 a { 2 . 1 A 2 2 \ MA 2 a M 2 
L = Vtt [-A 2 a(f) - —— [a 2 k + — + 4reV + 



2 V a 2 J y/a 2 + b 2 2^2b 

M 2 b 2 3H 2 A 2 a ( 1 A 2 2 \ 2 E 2 M 2 , , , 

+ — 7= + 4reV =-} , (17) 

8V2b 4 \a 2 J 2v/2& 3 

depending on the dynamical variables 0, re, A, a, b and their derivatives. The 
last two terms at ffT7|) contains the quantum corrections. 

Variation in <fi gives iV = 0, just reproducing the conservation of the high 
frequency quanta. Variation in re, the so-called chirp function, gives 

aa = 4rea 2 + 12# 2 re(l + 4reV) . (18) 

In contrast to the classical case where the chirp function is easily derived 
from a(t), in the quantum case re is given in terms of a as the solution of the 
third-degree equation ffTBl . Hence, expressing re in terms of a and a would 
give equations too cumbersome to be of any value. A reasonable alternative 
in the semiclassic limit is to solve ([TBI for re as a power series in H 2 . Other 
possibility is to regard (TIB"]) as a dynamical equation to be numerically solved 
for a(t). Both approaches will be considered in what follows. 
Combining independent variations in a and in A gives 

^y^-wfi^ +3 -^- imW ' (19) 

while varying b gives 

.. _ 1_ 2y/2 Ns z 3H 2 

S ~^~ M{a 2 + s A fl 2 + ~V ' ( ' 

where b = s 2 . 

Equations ( [1811201) form a complete system for the dynamical variables 
a, re and s and are the basis for the conclusions in the following. The system 
(H8H20p will be studied in the semiclassic limit, both for linear and nonlinear 
oscillations, and for arbitrary values of the quantum parameter H. 



3 Semiclassic oscillations 

It is of considerable interest to investigate the modifications in the Langmuir 
soliton induced by small quantum effects. When if is a small parameter, 
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(I18p can be solved approximately in powers of H 2 yielding 



d tt9 / 3d 3d 3 \ .„„ . 

k = H — - + , 21 

4a \4a 3 16a/ V ; 

disregarding higher-order corrections. Then one is left with a coupled, non- 
linear system of second-order equations for a and s, namely, (120]) and 

4 4Ma , /24 6d 2 12M 9Mad 2 \ 



a 3 (a 2 + s 4 ) 3 / 2 \a 5 a 3 a(a 2 + s 4 ) 3 / 2 (a 2 + s 4 ) 3 / 2 , 

Consistently, the energy (Tl6]) evaluated using (l2Tj) . 



W = l( iV 4 + a^ ) + ^ M2(i2 + ^ 
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is approximately constant, dTL/dt = 0(H 4 ) along trajectories of (|2Q]) - (|22]) . 
and can be used to check the accuracy of numerical schemes. 

It is relevant to check (r2Ul) - (l22l) for the linear stability of fixed points. For 
the Langmuir soliton fl9l fTU|) one have equal values for M and N. Since we 
are mainly interested in the role of quantum effects for the Langmuir soliton, 
in the rest of the section we set M = N. Since quantum effects are small, 
one can search for fixed points for the dynamical system as a power series in 
H 2 . An easy calculation then yields critical points at (a, s) = (a , Sq), with 

2y/2 3H 2 M 

a " = ir + ^r- < 24) 

/2V2\' /2 H 2 M 3 ' 2 , v 

*> = \-m) + ^J^' < 25 » 

disregarding 0(H 4 ) terms. In terms of the original variables a and b these 
fixed points corresponds to Gaussians of same width at the formal classical 
limit. Quantum corrections, however, introduce a disturbance: the width a 
of the Gaussian for the envelope electric field increases less than the width 
b associated to the density. Moreover, both characteristic lengths increase, 
pointing for a wave-packet spreading effect. 
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Considering small deviations ~ exp(iut) from the equilibrium point, one 
obtain 



a; 2 



± 



M 2 r i 

24 + 10M 2 - H 2 M 2 (18M 2 + 33) (26) 

128 ^ 

M 2 
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144 + 24M 2 + 25M 4 - H 2 M 2 {396 + 177M 2 + 90M 4 ) 



1/2 



For consistency one could also expand the square root at (1251) up to 0(H 2 ) 
terms, but this would result in a more cumbersome expression. 

Unstable linear oscillations corresponds to solutions with Im(uj) < 0. 
A straightforward algebra shows that such instabilities are impossible in the 
formal classical limit H = 0. In the quantum case, however, a careful analysis 
of ( J26l) shows that instabilities are possible when 

f(M 2 )= 144 + 24M 2 + 25M 4 

M J ~ M 2 (396 + 177M 2 + 90M 4 ) ' { ' 

This instability condition can be satisfied for small values of H. For instance, 
for M > 0.6, f(M 2 ) < 1 at flZZ}. Further increasing; M allows for smaller 
values of if. For M = 5, one has if > 1/10 for instability, an inequality 
which can be satisfied within the present semi classic context. In terms of 
the Langmuir soliton (19T TT0]) . M = 2y/2E , showing that large amplitude 
solitons are more influenced by quantum instabilities. Figure 1 shows the 
curve H 2 = f(M 2 ) separating stable- unstable regions. 

The existence of quantum instabilities not present for the classical system 
is a remarkable fact. It shows that the classical localized solution eventually 
is smeared out, since the width of the Gaussian is continuously increasing in 
time. This is a signature of wave-packet spreading and tunneling, because 
the classical Langmuir soliton is produced by particle trapping in the self- 
consistent electrostatic potential. Notice, however, that nonlinear effects can 
suppress the linear quantum instability. 

Unlike the classical case [26], the system (I20l) - (l22l) seems to be not de- 
scribed by a pseudo-potential function. This is due to the velocity depen- 
dence on the dynamical equation for a(t). Nevertheless, using the energy 
integral 7i one can get a rough estimate for an escape velocity. In the limit 
situation where the particle escapes, one has a = s = when a — > oo and 
s — > oo. From (1231) it gives Ji — 0. Now supposing an initial condition at 
the fixed point (I24H25P and taking d(0) = d = for simplicity, one finds a 



9 



1.5 



0.5 











2 

M 2 



Figure 1: The curve H 2 = f(M 2 ) 
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for the instability condition ( 1271) . Instabil 



escape velocity s(0) = s such that (for M = N) 

M 7H 2 M 3 



(28) 



u AV2 6AV2 

Once again, quantum effects act in a tunneling- like manner. Indeed, 
shows that a sufficiently large value of H can produce escaping of the par- 
ticle, no matter the value of the initial velocity Sq. The limiting value 
H = (4/v / 7)M" 1 for which s = can be achieved even for the semiclassic 
case for sufficiently high M. For instance, when M = 10, the particle even- 
tually escapes for H > 0.15, a moderate value. A nonzero value of do tends 
to produce even smaller values of H for escaping. 

Figures 2, 3 and 4 shows typical oscillations for the dynamical system 
p)]) - (l22]) . with M = N = 3, H = 0.3. The initial condition is at the fixed 
point and d(0) = 0. Also, so = 0.62. For such parameters, simulations shows 
unbounded motion for so = 0.64, which is much less than the classical escape 
velocity, Sq = 0.73, and in good agreement with the critical value 0.59 arising 
from the crude estimate ff28l . The Hamiltonian TC remains approximately 
constant at the value —0.10 along the run. Observe the different time-scales 
for a(t) and s(t). Taking a smaller value of Sq gives a more regular, quasi- 
periodic oscillation pattern, similar to the classical oscillations ■ However, 
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t 

Figure 2: Simulation for the semiclassic system (I20 |) -( l22|) showing a(t). Pa- 
rameters, M = N = 3, H = 0.3. Initial condition, (ao, so, a>o, so) = 
(1.52,1.36,0,0.62). 

notice that quantum effects leads to complicated trajectories even for equal 
values of the invariants M and N, approaching the critical value of So f° r 
unbounded motion (see figure 4). Several other runs shows that increasing 
the value of so increases the period and the amplitude of the oscillations. 
Direct comparison between the present simulations and those of the original 
quantum Zakharov system will be reported in a future work. 

4 Non perturbative full system 

When H is not a small parameter, one is not allowed to solve (|T8|) retaining 
only the first-order quantum correction. Therefore, one is left with the full 
system (I18l)- (l20l . for which some conclusions can be obtained. Even for the 
balanced case when M = N, it is not possible to get a closed-form solution 
for the fixed points, making difficult to derive general statements about lin- 
ear stability. However, simulations can be made for different values of H, 
starting from a fixed point numerically calculated. For M — N — 1, H — 5, 
an equilibrium is found for (k, a, s) = (0, 9.15, 3.53). Figure 5 shows a typical 
trajectory starting at this initial condition, with so = 0.2. Under the same 
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Figure 3: Simulation for the semiclassic system (I20 i) -( l22|) showing s(t). P 
rameters, M = N = 3, H = 0.3. Initial condition, (ao, so, ao, ^o) 
(1.52,1.36,0,0.62). 
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Figure 5: Simulation for the full dynamical system ( |T8l) -(l20l) showing a and 
s. Parameters, M = N = 1, H = 5. Initial condition at (k, a, s, s) = 
(0,9.15,3.53,0.20). 



parameters but with a smaller initial velocity produces quasi-periodic motion, 
as shown in figure 6, where so = 0.05. Similar simulations shows that for 
increasing H it becomes more difficult to get regular, quasi-periodic trajec- 
tories, pointing for instabilities of quantum nature. In addition, unbounded 
motion appears for smaller values of the initial velocity. 

Further results on the role of quantum effects can be obtained for the 
ultra-quantum case where we can neglect all terms at the right-hand sides of 
( TT8]) - (l20]) not containing H 2 . In such case where H 2 is large enough, one get 
the system 

a = 12# 2 k(1 + 4/tV)/a, (29) 
k = 3H 2 (l-lQK 4 a 8 )/a\ (30) 
s = 3H 2 /s 7 . (31) 

Equations ( 12911301) can be solved yielding 

2 _ a 3QH\t - tp) 2 9a 4 H\t-t ) 2 
a ~ %+ al ' K ~ (a* + 36H*(t - to) 2 ) 2 > {6Z) 

where ao and to are numerical constants. From ( |32|) the conclusion is that the 
width a of the envelope electric field tends to increase without bound in the 
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Figure 6: Simulation for the full dynamical system ( |T8l) -(l20l) showing a and 
s. Parameters, M = N = 1, H = 5. Initial condition at (k, a, s, s) = 
(0,9.15,3.53,0.05). 

ultra-quantum limit, while the chirp function k approaches zero. Similarly, 
inspection of fl3Tl) shows that the width b = s 2 of the density fluctuation 
increases without bound, since s = —dV/ds for a pseudo-potential V = 
H 2 /2s 6 having no bound-states. The results at the end of this section comes 
from the rough procedure of disregarding all terms not containing H 2 in 
the right-hand sides of (jT8]) - (!20l . They are consistent with the previous 
simulations, showing the destabilizing role of quantum effects. 

5 Conclusion 

The quantum Zakharov system was analyzed through a time-dependent Gaus- 
sian trial function method for an associated Lagrangian formalism. This 
extends to the quantum plasmas realm the results obtained for the classi- 
cal Zakharov system by a similar approach [25]. In contrast to the classical 
case, complicated trajectories and instabilities can be found even for the bal- 
anced case of equal values of the invariants M and N, corresponding to the 
low-frequency part of the mass and the number of high frequency quanta, 
respectively. Quantum effects plays a destabilizing role, yielding the ulti- 
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mate decaying of the Langmuir soliton whose properties are simulated by 
the time-dependent Gaussian ansatz. This is a signature of quantum effects 
(wave-packet spreading, tunneling), making more difficult the existence of 
coherent, localized solutions in quantum plasmas. Direct comparison be- 
tween the variational solutions of this work and numerical simulations of the 
original Zakharov system remains an open question. 
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